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Abstract 

The integrability (solvability via an associated single-valued linear 
problem) of a differential equation is closely related to the singular- 
ity structure of its solutions. In particular, there is strong evidence 
that all integrable equations have the Painleve property, that is, all 
solutions are single-valued around all movable singularities. In this 
expository article, we review methods for analysing such singularity 
structure. In particular, we describe well known techniques of non- 
linear regular-singular-type analysis, i.e. the Painleve tests for ordi- 
nary and partial differential equations. Then we discuss methods of 
obtaining sufficiency conditions for the Painleve property. Recently, 
extensions of irregular singularity analysis to nonlinear equations have 
been achieved. Also, new asymptotic limits of differential equations 
preserving the Painleve property have been found. We discuss these 
also. 
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1 Introduction 



A differential equation is said to be integrable if it is solvable (for a sufficiently 
large class of initial data) via an associated (single-valued) linear problem. 
A famous example is the Korteweg-de Vries equation (KdV), 

Ut + QUU-^ + Uxxx = 0, (1.1) 

where the subscripts denote partial differentiation. 

The KdV equation was discovered to be integrable by Gardner, Greene, 
Kruskal, and Miura |21|. (Its method of solution is called the inverse scatter- 
ing transform (1ST) method; see the paper by Mark Ablowitz in the present 
collection.) Since this discovery, a large collection of nonlinear equations 
(see |]ll) has been identified to be integrable. These range over many di- 
mensions and include not just partial differential equations (PDEs) but also 
differential-difference equations, integro-differential equations, and ordinary 
differential equations (ODEs). 

Six classical nonlinear second-order ODEs called the Painleve equations 
are prototypical examples of integrable ODEs. They possess a characteristic 
singularity structure i.e. all movable singularities of all solutions are poles. 
Movable here means that the singularity's position varies as a function of 
initial values. A differential equation is said to have the Painleve property 
if all solutions are single-valued around all movable singularities. (See com- 
ments below and in Section 2 on variations of this definition.) Thefore, the 
Painleve equations possess the Painleve property. Painleve [Q, Gambler 
pO| , and R. Fuchs ||l9l identified these equations (under some mild condi- 
tions) as the only ones (of second order and first degree) with the Painleve 
property whose general solutions are new transcendental functions. 

Integrable equations are rare. Perturbation of such equations generally 
destroys their integrability. On the other hand, any constructive method of 
identifying the integrability of a given system contains severe shortcomings. 
The problem is that if a suitable associated linear problem cannot be found it 
is unclear whether the fault lies with the lack of integrability of the nonlinear 
system or with the lack of ingenuity of the investigator. So the identification 
of integrability has come to rely on other evidence, such as numerical studies 
and the singularity structure of the system. 

There is strong evidence 0, ^ that the integrability of a nonlinear sytem 
is intimately related to the singularity structure admitted by the system in its 
solutions. Dense multi-valuedness (branching) around movable singularities 
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of solutions is an indicator of nonintegrability|p^ . The Painleve property 
excludes such branching and has been proposed as a pointer to integrability. 

The complex singularity structure of solutions was first used by Kowalevs- 
kaya ^ to identify an integrable case of the equations of motion for a 
rotating top. Eighty eight years later, this connection was reobserved in 
the context of integrable PDEs by Ablowitz and Segur 0, and Ablowitz, 
Ramani, and Segur Their observations led to the following conjecture. 

The ARS Conjecture: Any ODE which arises as a reduc- 
tion of an integrable PDF possesses the Painleve property, possi- 
bly after a transformation of variables. 

For example, the sine-Gordon equation 

Ma;t = sinM, (1.2) 

which is well known to be integrable admits the simple scaling sym- 

metry 

X Ax, t ^ A^^t. (1.3) 

To find a reduction with respect to the symmetry ( |1.3| ), restrict to the sub- 
space of solutions that is invariant under ( p..3[ ) by introducing new variables 
z, w such that 

u{x,t) = w{z), z = xt. 

This gives 

zw" + w' = sin If, (1.4) 

where the prime denotes differentiation with respect to z. To investigate 
the Painleve property, this equation must first be transformed to one that 
is rational (or possibly algebraic) in w. (Otherwise, the nonlinear analogue 
of Frobenius analysis used to investigate the Painleve property cannot find 
a leading-order term to get started. See Section 2.) Introduce the new 
dependent variable y := exp(iw). Then equation ( |1.4| ) becomes 

z{yy" - y'^) + yy' = ^y{y^ - l). 

This equation (a special case of the third Painleve equation) can be shown 
to have the Painleve property. (See Section 2.) 

There is now an overwhelming body of evidence for the ARS conjecture. A 
version that is directly applicable to PDEs, rather than their reductions, was 
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given by Weiss, Tabor, and Carnevale (WTC)|^. The ARS conjecture and 
its variant by WTC are now taken to be almost self-evident because they have 
been formally verified for every known analytic soliton equation |5^, |51[] 
(where analytic means that the equation is, or may be written to be, locally 
analytic in the dependent variable and its derivatives). Previously unknown 
integrable versions of the sohton equations |^ |15|, ^ have been identified 
by the use of the conjecture. The conjecture has also been extrapolated to 
identify integrable ODEs p7[] . 

Rigorous results supporting the conjecture exist for ODEs with special 
symmetries or symplectic structure pO, E^. Necessary conditions for pos- 



sessing the Painleve property |17, 29 1 have also been derived for general 



semilinear analytic second-order PDEs. No new integrable PDEs were found 
— suggesting that in this class at least, there is a one-to-one correspondence 
(modulo allowable transformations) between integrable equations and those 
possessing the Painleve property. Moreover, proofs of weakened versions of 
the conjecture exist |^ These results point strongly to the truth of the 
ARS conjecture. Nevertheless, the conjecture has not yet been proved. 

The main aim of this paper is to describe methods for investigating the 
singularity structure of solutions of ODEs and PDEs. These may be divided 
into two classes, those that parallel methods for analysing regular singular 
points and those that parallel techniques for irregular singular points of linear 
ODEs. 

The first class of methods has been widely used formally. The most 
popular procedure is to expand every solution of the differential equation of 
interest in an infinite series near a movable singularity of the equation |^ , 
i.e. the solution u{z) is expanded as 

oo 

uiz) = Y.a^iz-zor^'' (1.5) 

n=0 

where zq is the arbitrary location of a singularity and p is the leading power 
that needs to be found. Such an expansion is often called a Painleve expan- 
sion. The equation is assumed to have the Painleve property if the series 
is self-consistent, single- valued, and contains a sufficient number of degrees 
of freedom to describe all possible solutions or the general solution. These 
demands yield necessary conditions for the Painleve property to hold. The 
series and the expansion techniques are analogues of the usual Frobenius 
(or Fuchsian) expansion procedure for linear ODEs. This procedure was 
extended to PDEs by 
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These techniques are, in general, not sufficient to prove that a differential 
equation has the Painleve property. For example, even if the only possible 
formal solutions are Laurent series, the poles indicated by these series may 
accumulate elsewhere to give rise to a worse (branched) singularity. 

Painleve gave sufficient conditions to show that his eponymous equations 
have the Painleve property. However, his proof is not widely understood. 
We describe briefly here an alternative, direct, method of proof due to Joshi 



and Kruskal |Q. To gain sufficiency, we showed that the solutions of the 
Painleve equations possess convergent Laurent expansions around every mov- 
able singularity, and moreover, (in any given bounded region) the radius of 
convergence of each series is uniformly bounded below. In other words, the 
poles of any solution cannot coalesce to form a more complicated singularity 
elsewhere (in the finite plane). 

For nonlinear PDEs, the question of how to get sufficient conditions for 
the Painleve property is still open. Nevertheless, partial results are now 
known. These make the WTC analogue of the Frobenius method rigorous 
and go some way toward proving that a given PDE has the Painleve property. 



We describe the results due to Joshi and Petersen [051 361 and Joshi and 



Srinivasan |^8[ in section 2. An alternative approach to the convergence of 
the Painleve expansions for PDEs has also been developed by Kichenassamy 
and Littman 0, |T|]. 

Another difficulty with the analysis of singularity structure is that the 
Painleve expansions can miss some solutions. This may happen, for exam- 
ple, when the number of degrees of freedom in the series is less than the 
order of the differential equation. Perturbations of such series often reveal 
that the missing degrees of freedom lie in terms that occur (paradoxically) 
before the leading term. For reasons explained in Section 2, such terms are 
called negative resonances. In other cases, perturbations reveal no additional 
degrees of freedom at all. We call the latter series defective. 

How can we deduce the singular behaviours of solutions that are missed by 
the Painleve expansions? We provide an answer based on irregular-singular 
analysis for linear ODEs and illustrate it through two important examples. 
The ffist example is the Chazy equation |l[ , a third-order ODE, whose gen- 
eral solutions have movable natural barriers. The Painleve expansion of the 
solution of the Chazy equation contains only two arbitrary constants. The 
second example is a fourth-order ODE ffist studied by Bureau. This example 
has two families of Painleve expansions, one of which has negative resonances 
and one that is defective. In section 3, we show how the Painleve expansions 
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can be extended through exponential (or WKB-type) perturbations. 

Conte, Fordy, and Pickering |jl6[ have followed an alternative approach. 



Their perturbations of Painleve expansions involve Laurent series with no 
leading term (i.e. an infinite number of negative powers). As pointed out 
by one of us, this is well defined only in an annulus where the expansion 
variable is lower-bounded away from the singularity. Conte et al. demand 
that each term of such a perturbation must be single-valued. Therefore, their 
procedure requires a possibly infinite number of conditions to be checked for 
the Painleve property. Our approach overcomes this problem. 

For linear differential equations, in general, the analysis near an irregular 
singularity yields asymptotic results, i.e. asymptotic behaviours along with 
their domains of asymptotic validity near the singularity. The latter is crucial 
in this description. For example, it is well known that the Airy function Ai{x) 
which solves the ODE 

y" = xy, 

has the asymptotic behaviour 

Ai{x) ~ , exp(— 2a:^^^/3) as |x| oo, | argx| < vr 

near the irregular singular point at infinity. (See p, 0.) Note that the 
asymptotic behaviour of Ai{x) is apparently multivalued but the function 
itself is single-valued everywhere. (In fact, Ai{x) is entire, i.e. it is analytic 
throughout the whole complex x-plane.) 

The resolution of this apparent paradox lies in the angular width of the 
sector of validity of the above behaviour, which is strictly less than 27r. To 
describe Ai{x) in the whole plane, we need its asymptotic behaviour in re- 
gions that include the line | arg(x)| = vr. (Such behaviours are well known. 
See e.g. These, together with the behaviour given above, show that the 
analytic continuation of Ai{x) along a large closed curve around infinity is 
single-valued. Therefore, the global asymptotic description is not actually 
multivalued. 

On the other hand, suppose an asymptotic behaviour is multivalued and 
its sector of validity extends further than 2n. Then there are (at least) two 
asymptotic descriptions of a solution at the same place (near an irregular 
singularity). This violates the uniqueness of asymptotic description of a so- 
lution, unless the solution is itself multivalued. Therefore, such an asymptotic 
behaviour is an indication that the solution cannot satisfy the requirements 
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of the Painleve property. The Bureau equation we study in Section 3 provides 
an example of this case. 

Such results form an important extension of the usual tests for the Painleve 
property. However, there is no denying that many fundamental questions re- 
main open in this area, even at a formal level. For example, the Painleve 
property is easily destroyed by straightforward transformations of the de- 
pendent variable(s). (E.g. A solution u{z) of an ODE with movable simple 
poles is transformed to a function w{z) with movable branch points under 
u I— > w"^.) An extension of the Painleve property called the poly-Painleve 



property has been proposed by Kruskal p5|, p8| to overcome these difficul- 



ties. It allows solutions to be branched around movable singularities so long 
as a solution is not densely valued at a point. However, such developments 
lie outside the scope of this paper and we refer the reader to |^ for further 
details and references. 

Other major problems remain. One is to extend the classification work 
of Painleve and his colleagues to other classes of differential equations. Cos- 
grove has accomplished the most comprehensive extensions in recent times 



T7| . The universal method of classification called the a-method is based on 



asymptotic ideas (see Section 2). Asymptotic limits of differential equations 
can illuminate such studies. 

The Painleve equations are well known to have asymptotic limits to other 
equations with the Painleve property. These limits are called coalescence 
limits because singularities of the equation merge under the limit. In Section 
4, we describe the well known coalescence limits of the Painleve equations 
and show that these limits also occur for integrable PDEs. 

Throughout this article, solutions of differential equations assumed to be 
complex-valued functions of complex variables. 



2 Nonlinear- Regular- Singular Analysis 

In this section, we survey the main techniques used to study the Painleve 
property. These range from the a-method to the widely used formal test 
known as the Painleve test. 

Consider the second-order linear ODE 

u"{z) +p{z)u{z) + q{z)u{z) = 0, 

where primes denote differentiation with respect to z. Fuchs' theorem 
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states that u can only be singular (nonanalytic) at points where p and q 
arc singular. Such singularities are called fixed because their positions are 
determined d priori (before solving the equation) and their locations remain 
unchanged throughout the space of all possible solutions. 

However, fixed singularities are not the only possibilities for nonlinear 
equations. Consider the Riccati equation 



where Zq is an arbitrary constant. If, for example, the initial condition is 
u{0) — 1, then zq — —1. If the initial condition is changed to u{0) = 2, then 
zq moves to zo — —1/2. In other words, the location of the singularity at zo 
moves with initial conditions. Such singularities are called movable. 

Nonlinear equations exhibit a vast range of types of movable singularities. 
Some examples are given in the Table 1 (where k and Zq are arbitrary constant 
parameters) . 



u"{z) + u'^{z) = 0. 



It has the general solution 
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Table 1: Examples of Possible Singular Behaviour 




Equation 


General Solution 


Singularity Type 


1. 


y' + y^ = 


y={z- zo)-'^ 


simple pole 


2. 


2yy' = 1 


y^y/z-Zo 


branch point 


3. 


y" + 7/'2 = 


y = ln{z - Zo) + k 


logarithmic 
branch point 


4. 


vv" 

+ y'\y/y'-l^ 
= 


y ^kexp{[z - zo]~'^) 


isolated 
essential 
singularity 


5. 


(1 + y')y" 

= 


y = tan (ln(A;[2; - ^o])) 


nonisolated 

essential 
singularity 


6. 


{y" + y'y'f 

= i/V'(%' + 2 


y — A;tan[A;^(a; — xq)] 
or 

y=((4/3)/(x-Xo))'^' 


pole 

branch point 



A normahzed ODE, i.e. one that is solved for the highest derivative, such 

as 

y^-)^F{y^^-'\...,y',y,z), (2.6) 

gives rise to possible singularities in its solutions wherever F becomes sin- 
gular. (Where F is analytic, and regular initial data are given, standard 
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theorems show that an analytic solution must exist.) Note that these sin- 
gularities may include the points at infinity in y (or its derivatives) and 
which we denote by y = oo (or = oo etc), z = oo. The singularities of 
F, therefore, denote possible singularities in the solution(s). They may be 
divided into two classes: those given by values of z alone and those involving 
values of y or its derivatives. The former are determined d priori for all 
solutions. Therefore, they can only give rise to fixed singularities. To find 
movable singularities, we therefore need to investigate the singular values of 
F that involve y (or its derivatives). Similar statements can be made in the 
case of PDEs. 

For example, the Riccati equation 



has a right side F with two singularities given by 2; = and y = 00. The 
general solution is 



It is clear that z = is a fixed singularity (it stays the same for all ini- 
tial conditions) whereas zq denotes a movable singularity where y becomes 
unbounded. 

Singularities of nonlinear ODEs need not only occur at points where y 
is unbounded. Example 2 of Table 1 indicates possible movable singularities 
at points where y — 0. The solution shows that these are actually movable 
branch points. 

These considerations show that singular values of the normalized differen- 
tial equation lie at the base of the solutions' singularity structure. Techniques 
for investigating singularity structure usually focus on these singular values. 

In the first three subsections below, we describe common definitions of 
the Painlcvc property, and the two major techniques known as the a-method 
and the Painleve test for deriving necessary conditions of the Painleve prop- 
erty. In the subsequent three subsections below, we discuss the need for 
sufficiency conditions, the direct method of proving the Painleve property, 
and convergence-type results for PDEs. 

2.1 The Painleve Property 

The actual definition of the Painleve property has been subject to some 
variation. There are three definitions in the literature. 



y' = -y^/z =: F{y,z). 
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Definition 2.1 An ODE is said to possess 

1. the specialized Painleve property if all movable singularities of all solu- 
tions are poles. 

2. the Painleve property if all solutions are single-valued around all movable 

singularities. 

3. the generalized Painleve property if the general solution is single-valued 
around all movable singularities. 

(The qualifiers "specialized" and "generalized" are not usually used in the 
hterature.) The first property defined above clearly imphes the others. This 
property was also the first one investigated (by ARS) in recent times. It is 
the property possessed by the six Painleve equations. 

The second, more general, definition above is the one used by Painleve 
in his work on the classification of ODEs. It allows, for example, movable 
unbranched essential singularities in any solution. Of the examples in Table 
1, equations 1 and 4 have the Painleve property; equation 1 also has the spe- 
cialized Painleve property. The remaining equations have neither property. 

The third property is the most recently proposed variation, although 
there is evidence that Chazy assumed it in investigating ODEs of higher 
(> 1) degree or order (> 2). The sixth example given in Table 1 satisfies 
neither of the first two properties above because the special solution (4/3/ {x— 
xq))^^^ has movable branch points around which the solution is multivalued. 
However, it does satisfy the generalized Painleve property because the general 
solution ki3xi[k^{x — Xq)] is meromorphic. 

Most of the known techniques for investigating the Painleve property 
have their origin in the classical work of Painleve and his colleagues. They 
classified ODEs of the form 

u" ^ F{z]u,u'), (2.7) 

where F is rational in u and u' and analytic in according to whether or 
not they possess the Painleve property. 

They discovered that every equation possessing the Painleve property 
could either be solved in terms of known functions (trigonometric functions, 
elliptic functions, solutions of linear ODEs, etc.) or transformed to one of 
the six equations now called the Painleve equations {Pj-Pvi)- They have 
standard forms that are listed below. (They are representatives of equivalence 
classes under mobius transformations.) Their general solutions are higher 
transcendental functions. 



11 



The Painleve Equations 



u" = 6u^ + z, 

u" = 2u^ + zu + a, 

u" = -M'^--u' + -(aM^ + /3)+7M^ + -, 
u z z u 

u" = —u'^ + -u^ + Azu^ + 2(z'^-a)u+-, 
2u 2 u 

u = — + -\u --U+ — « + - + — + — , 

I2u U—1) z z'^ \ u z u — \ 



2 Vu u — 1 u — z ) Iz u — 1 u — z 

u{u-l){u-z) \ (3z 7(z - 1) 6z{z - 1) ] 

Z2(^_l)2 {U-If {U- Zf ]' 

Two main procedures were used in tliis work. Tlie first is known as the 
a-method and the second is now called Painleve analysis. Painleve described 
the a-method in the following way. 

Considerons un equation differentielle dont le coefficient differentiel 
est une fonction (holomorphe pour a = 0) d'un parametre a. Si 
I'equation a ses points critiques fixes pour a quelconque (mais 
0), il en est de meme, a fortiori pour a = 0, et le developpement 
de I'integrale y{x), suivant les puissances de a, a comme coefficients 
des fonctions de x a points critiques fixes. 



(This extract is taken from footnote 3 on p. 11 of ||55[. In Painleve's terminol- 
ogy, a critical point of a solution is a point around which it is multivalued.) 
In other words, suppose a parameter a can be introduced into an ODE in 
such a way that it is analytic for a = 0. Then if the ODE has the Painleve 
property for a 7^ 0, it must also have this property for a = 0. We illustrate 
this method for the classification problem for first-order ODEs below. 

The second procedure, called Painleve analysis, is a method of examining 
the solution through formal expansions in neighbourhoods of singularities of 
the ODE. In particular, the procedure focusses on formal series expansions 
of the solution(s) in neighbourhoods of generic (arbitrary) points (not equal 
to fixed singularities). The series expansion is based on Frobenius analysis 
and usually takes the form given by Eqn (|1.5| ). 
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As mentioned in the Introduction, this procedure was extended to PDEs 
by WTC (Weiss, Tabor, and Carnevale) [^. For PDEs, the above defini- 
tions of the Painleve property continue to hold under the interpretation that 
a movable singularity means noncharacteristic analytic movable singularity 
manifold. 

A noncharacteristic manifold for a given PDE is a surface on which we 
can freely specify Cauchy data. The linear wave equation, 

utt - u^x = 0, (2.8) 

has the general solution 

u{x, t) = fit — x) + git + x), 

where / and g are arbitrary. By a suitable choice of / and g we can construct 
a solution u with any type of singularity on the curves t — x = ki, t + x = k2, 
for arbitrary constants ki, and /c2- These lines are characteristic manifolds 
for equation ( p.8|) . This example illustrates why the Painleve property says 
nothing about the singular behaviour of solutions on characteristic singularity 
manifolds. 

The WTC procedure is to expand the solutions uix,t) of a PDE as 

oo 

uix,t) = J2unix,t)^^+P, (2.9) 

n=0 

near a noncharacteristic analytic movable singularity manifold given by $ = 
0. (This extends in the obvious way for functions of more than two variables.) 
The actual expansion can be simplified by using information specific to the 
PDE about its characteristic directions. For example, for the KdV equation, 
noncharacteristic means that $2; 7^ 0. Hence by using the implicit function 
theorem near the singularity manifold, we can write 

<!>{x,t) =x-at), 

where ^{t) is an arbitrary function. This is explored further in subsection 
2.3.2 below. 

In some cases, the form of the series Eqn( |1.5D (or (( p.9|) ) needs modifica- 
tion. A simple example of this is the ODE 

u"' = 2{u'f + l. (2.10) 
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Here v = u' is a Jacobian elliptic function with simple poles of residue ±1. 
(See 0.) Hence a series expansion of u{z) around such a singularity Zq, say, 
must start with ±log(2; — zo)- The remainder of the series is a power series 
expansion in powers oi z — zq. In such cases, the Painleve property holds for 
the new variable v. 



2.2 The a-Method 

In this section, we illustrate the a-method by using it to find all ODEs of 
the form 

, P{z, u) 

«' = 75rM 2.11 

possessing the Painleve property, where P and Q are analytic in z and poly- 
nomial in u (with no common factors). The first step of the a- method is to 
introduce a small parameter a through a transformation of variables in such 
a way that the resulting ODE is analytic in a. However, the transformation 
must be suitably chosen so that the limit a — > allows us to focus on a 
movable singularity. This is crucial for deducing necessary conditions for the 
Painleve property. 

We accomplish this by using dominant balances of the ODE near such 
a singularity. (See for a definition and discussion of the method of 

dominant balances.) 

If Q has a zero of multiplicity m aX u = a{z) then, after performing the 
transformation u{z) ^ u{z) + a{z), equation ( p.ll| ) has the form 

u'^u' = f{z,u), (2.12) 

where / is analytic in m at m = and /(-2,0) ^ (since P and Q have 
no common factors). Choose zq so that n := /(zo,0) ^ and define the 
transformation 

u{z) = aU{Z), z = Zq + a'^Z, 

where n is yet to be determined and a is a small (but nonzero) parameter. 
Note that this is designed to focus on solutions that become close to the 
singular value m = of the equation somewhere in the z-plane. 



Equation ( |2.12| ) then becomes 



^m+i-n^m^ = f {zq + a^'Z, aU) = K + 0{a). 
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This equation has a dominant balance when n = m + 1. In this case the hmit 
as a ^ gives 

which has the exact solution 

U{Z) = {{m + 1)kZ + C}^ , 

where C is a constant of integration. This solution has a movable branch 
point at Z = —C/{K,{m + l)) for all m > 0. Therefore, Eqn cannot 
possess the Painleve property unless m = 0, i.e. Q must be independent of 
u. That is, to possess the Painleve property, Eqn must necessarily be 

of the form 

u' = ao{z) + ai{z)u + a2{z)u^ + ■ ■ ■ + a]\f{z)u^ , (2-13) 

for some nonnegative integer N. 

The standard theorems of existence and uniqueness fail for this equation 
wherever u becomes unbounded. To investigate what happens in this case, 
we transform to v := 1/u. (Note that the Painleve property is invariant 
under such a transformation.) Eqn (|2.13|) then becomes 

v' + ao{z)v'^ + ai{z)v + a2{z) + a3{z)v^^ + ■ ■ ■ + a]^{z)v'^~'^ = 0. 

But this equation is of the form and so it can only possess the Painleve 

property if = 2. 

In summary, for Eqn ( |2.11| ) to possess the Painleve property, it must 
necessarily be a Riccati equation: 

u' = ao{z) + ai{z)u + a2{z)u'^ . (2-14) 

To show that this is also sufficient, consider the transformation 

1 w' 



u 



a2{z) w 
which linearizes Eqn ( |2.14|) 

a2w" — + aia2)w' + a^alw = 0. 

By Fuchs theorem the singularities of any solution w can only occur 
at the singularities of (a'g + aia2)/a2 or aoa2- These are fixed singularities. 
Hence the only movable singularities of u occur at the zeroes of w. Since w 
is analytic at its zeroes, it follows that u is meromorphic around such points. 
That is, Eqn ( 2.14|) has the Painleve property. 
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2.3 The Painleve Test 



Here we illustrate the widely used formal tests for the Painleve property for 
ODEs and PDEs by using examples. 

2.3.1 ODEs 

Consider a class of ODEs given by 

m" = 6u" + /(2), (2.15) 

where / is (locally) analytic and n > 1 is an integer (the cases n = or 1 
correspond to linear equations). 

Standard theorems that yield analytic solutions fail for this equation 
wherever the right side becomes singular, i.e. where either f{z) or u be- 
comes unbounded. We concentrate on the second possibility to find movable 
singularities. This means that the hypothesized expansion Eqn( |1.5D must 
start with a term that blows up at zq. To find this term, substitute 

U{z) ~ Co{z - Z^Y, Z Zq, 

where 3?(p) < 0, Cq 7^ 0, into Eqn (|2.15|) . This gives the dominant equation 



cop(p -i){z- z^y-^ + o ((^ - z^y-^) = Qcliz - z^r + o[{z- z,r'p+\ 

(2.16) 

The largest terms here must balance each other (otherwise there is no such 
solution). Since Cq 7^ and p 7^ or 1, we get 

_9 

p — 2 = np, =^ p 



n 



If p is not an integer, then u is branched at zq. Hence, the only n > 1 for 
which equation (p.l5|) can possess the Painleve property are n = 2 and n = 3. 



We will only consider the case n = 2 here for conciseness. The case n = 3 
is similar. (The reader is urged to retrace the following steps for the case 
n = 3.) 

If n = 2 then p = —2 (which is consistent with our assumption that 
3fJ(p) < 0). Then equation ( |2.16D becomes 



6co{z-zo)-'' = 6cl + 0({z-zo^-' 



5 
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which gives Cq = 1. Hence the hypothesized series expansion for u has the 
form 



ulz 



£ c„(z - zo)"-'. (2.17) 



n=0 



The function f{z) can also be expanded in a power series in z — zq because 
by assumption it is analytic. Doing so and substituting expansion (|2.17| ) into 
eqn (|2.15|) gives 

oo 

^(^_2)(n-3)c„(z-Zo)"-' 

n=0 

oo oo -| 

= 6 y: c.fi,{z - z,y-^^-' + y: —f^"'\^o)iz - z,r 

i,j=0 m=0 

= 6iz - zq)-^ + Uciiz - zo)~^ 

+Q{cl + 2c2){z - zo)-^ + 12(c3 + ciC2){z - zo)-^ 

oo f n 1 

+ E 6 E cn^c._ + 7— :7Ty/^"-'H^o) {z - z,r-\ 
n=4 I m=o yn ) 

Equating coefficients of like powers of {z — Zq) we get Ci = 0, C2 = 0, C3 = 0, 
and 

-2){n- 3)Cn = 6 E CmCn-m + , \ , , /^""'H-go) . (^ > 4). 

(n-4)! 

Note that c„ appears on both sides of this equation. Solving for c„, we find 

n— 1 Y 

{n + \){n - 6)c„ = 6 E ^^mCn-m + 7 7TT/^"~'n^o). (2.18) 

^1 (n-4)! 

For each n 7^ 6, this relation defines c„ in terms of {cm}o<m<n- However, 
for n = 6, the coefficient of c„ vanishes and equation ( p.l8| ) fails to define 
Cg. If the right side also vanishes, Cg is arbitrary. However, if the right side 
does not vanish there is a contradiction which implies that the series ( p.l7| ) 
cannot be a formal solution of eqn( p.l5| ). 

In that second case, the expansion can be modified to yield a formal 
solution by inserting appropriate logarithmic terms starting at the index n = 
6. (This is also the case for Frobenius expansions when the indicial exponents 
differ by an integer - see IQ. See |5l[] also for a rigorous study of equations 
admitting such algebraico- logarithmic expansions in several variables.) In 
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such a case, logarithmic terms appear infinitely often in the expansion and 
cannot be transformed away (as in Eqn( |2.1(i| ) above). They therefore indicate 
multivaluedness around movable singularities. 

That is, Eqn( p.l5D fails the Painleve test unless the right side of Eqn( p.l8|) 



vanishes at n = 6. This condition reduces to 

f'izo) = 0. 

However, since zq is arbitrary, this implies /" = 0. That is, f{z) = az + b 
for some constants a, b. If a = 0, this equation can be solved in terms of 
(Weierstrass) elliptic functions. Otherwise, translating z and rescaling u and 
z gives 

u" = 6u^ + z, (2.19) 

which is the first Painleve equation. 

The index of the free coefficient, cq, in the above expansion is called a 
resonance. The expansion contains two arbitrary constants, Cq and zq, which 
indicates that it captures the generic singular behaviour of a solution (because 
the equation is second order). 

There is a standard method for finding the location of resonances which 
avoids calculation of all previous coefficients. We illustrate this method here 
for Pj. After determining the leading order behaviour, substitute the per- 
turbation 

u^{z- zo)-^ + --- + I3{z- zoY''^ 

where r > into equation (|2.19| ). Here f3 plays the role of the arbitrary 
coefficient. To find a resonance r, we collect terms in the equation that are 
linear in f3 and demand that the coefficient of (3 vanishes. This is equivalent 
to demanding that /3 be free. The resulting equation 

(r + l)(r-6) = 0, 

is called the resonance equation and is precisely the coefficient of on the 
left side of equation ( ^.18| ). 

The positive root r = 6 is precisely the resonance we found earlier. The 
negative root r = — 1, often called the universal resonance, corresponds to 
the translation freedom in zq. (Consider zq Zq + e. Taking |e| < \z — zq\, 
and expanding in e shows that r = — 1 does correspond to an arbitrary 
perturbation.) 
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Note, however, that r = — 1 is not always a resonance. For example, 
consider an expansion that starts with a nonzero constant term, such as 

1 + ai{z - zq) + 

Perturbation of zq does not add a term corresponding to a simple pole to 
this expansion. 

If any resonance is not an integer, then the equation fails the Painleve 
test. The role played by other negative integer resonances is not fully under- 
stood. We explore this issue further through irregular singular point theory 
in Section 3. 

For each resonance, the resonance condition needs to be verified, i.e. 
that the equation at that index is consistent. These give rise to necessary 
conditions for the Painleve property. If all nonnegative resonance conditions 
are satisfied and all formal solutions around all generic arbitrary points zq 
are meromorphic, the equation is said to pass the Painleve test. 

This procedure needs to be carried out for evey possible singularity of the 
normalized equation. For example, the sixth Painleve equation, Pyi, has four 
singular values in u, i.e. u = 0,1, z and oo. The expansion procedure outlined 
above needs to be carried out around arbitrary points where u approaches 
each such singular value. (Table 2 in Section 4 lists all singular values of the 
Painleve equations.) 

2.3.2 PDEs 

In this subsection, we illustrate the WTC series expansion technique with an 
example. Consider the variable coefficient KdV equation, 

ut + f(t)uu^ + g(t)u^^^ = 0. (2.20) 

Let cj){x,t) be an arbitrary holomorphic function such that S := {{x,t) : 
0(x,t)} = is noncharcteristic. The fact that 5* is noncharacteristic for 
equation ( |2.20| ) means that 

0. ^ (2.21) 

on S. By the implicit function theorem, we have locally 4'{x,t) = x — ^{t) 
for some arbitrary function ^{t). 

We begin by substituting an expansion of the form 

oo 

u{x,t) = Y.Unit)r^'' (2.22) 

n=0 



19 



into equation ( 2.20 ). The leading order terms give a 
cients of powers of gives 



-2. Equating coeffi- 



= 


uo = -I2g/f, 


(2.23) 


n = \ 


ui = 0, 


(2.24) 


n = 2 


U2 = e/f 


(2.25) 


n = 3 


U3 = Mo/(/Uo), 


(2.26) 


n> 4 


{n + l){n — 4:){n — 6)gUn 


(2.27) 



n— 4 

■/ + l)M„-fc-3Mfc+3 

fc=0 

+ {n - A)CUn-2 - <_3 



(2.28) 
(2.29) 



Arbitrary coefficients can enter at n = 4, 6 if the recursion relation is consis- 
tent. Consistency at n = 6 is equivalent to 

2 

I / ?;■„ \ 

0. 



Ur 



1 



g{t) = f{t)\ao / f{s)ds + bo 



This implies that 



where oq and bo are arbitrary constants. In this case, Eqn( |2.20D can be 
transformed exactly to the KdV equation (see Grimshaw [p^). 

In particular, for the usual form of the KdV (/(t) = 6, g(t) = 1) we have 
the formal series expansion 

r(t) 



u{x, t) 



+ 



+ M4(t)r 



/,2 ■ g 



(2.30) 



Questions of convergence and well-posedness, i.e. continuity of the solution 
as the arbitrary functions (^(t), U4{t), UQ{t)) vary, are discussed in Section 
2.6 below. 



2.4 Necessary versus sufficient conditions for the Painleve 
property 

The methods we described above can only yield necessary conditions for the 
Painleve property. Here we illustrate this point with an example (due to 
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Painleve) which does not possess the Painleve property, but for which the 
Painleve test indicates only meromorphic solutions. 
Consider the ODE 



1 + u^)u + (1 - 2u)u'^ = (2.31) 



(see Ince The singularities of this equation are u = ±z, u = oo and 

u' = oo. Series expansions can be developed for solutions exhibiting each 
of the above singular behaviours and the equation passes the Painleve test. 
This equation, however, has the general solution 

u{z) = tan{log[A;(2; — zq)]}, 

where k and zq are constants. For k 0, u has poles at 

z = Zq + k~^ exp{ — (n + l/2)7r} 

for every integer n. These poles accumulate at the movable point Zq, giv- 
ing rise to a movable branched nonisolated essential singularity there. This 
example clearly illustrates the fact that passing the Painleve test is not a 
guarantee that the equation actually possesses the Painleve property. 
This danger arises also for PDEs. The PDE 

Wt = {1+ W^)Wxx + (1 - 2w)Wa;, 

under the assumption 

W7(x, t) =: u{x) 

reduces to the ODE above. 

To be certain that a given differential equation possesses the Painleve 
property, we must either solve it explicitly or implicitly (possibly through 
transformations to other equations known to have the property), or develop 
tests for sufficiency. Most results in the literature rely on the former ap- 
proach. In the next section, we develop a method for testing sufficiency. 

2.5 A Direct Proof of the Painleve Property for ODEs 

In this section we outline a direct proof given in (Joshi and Kruskal ||3l|, 
p4[ ) that the Painleve equations indeed possess the Painleve property. The 
proof is based on the well known Picard iteration method (used to prove 
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the standard theorems of existence and analyticity of solutions near regular 
points) modified to apply near singular points of the Painleve equations. A 
recommended simple example for understanding the method of proof is 



u 



6u^ + l 



which is solved by Weierstrass elliptic functions. 

Consider the initial value problem for each of the six Painleve equations 
with bounded data for u and u' given at an ordinary point Zi G C (The 
point zi cannot equal for Pm — Pyi, 1 for Py or Pyi, or Zi for Pyi — see 
Table 2.) 



Table 2: Fixed and Movable Singularities 




of the Painleve Equations 


Equation 


Fixed Singularities 


Movable Singularities 




(2;- value) 


(w-value) 


Pi 


{00} 


{00} 


Pii 


{00} 


{00} 


Piii 


{0,00} 


{0,00} 


Piv 


{00} 


{0,00} 


Pv 


{0,00} 


{0,1,00} 


Pvi 


{0,1,00} 


{0,1, 2;, 00} 



Standard theorems yield a (unique) solution U in any region in which the 
Lipschitz condition holds. However, they fail where the right side becomes 
unbounded, i.e. at its singular values. (See e.g. |]14|-) Since our purpose 
is to study the behaviour of the solution near its movable singularities, and 
these lie in the finite plane, we confine our attention to an arbitrarily large 
but bounded disk \z\ < B (where B is real and say > 1). For Pm, Py, and 
Pvi this must be punctured at the finite fixed singularities. Henceforth we 
concentrate on Pj for simplicity. 

The ball \z\ < B contains two types of regions. Around each movable 
singularity, we select a neighbourhood where the largest terms in the equation 
are sufficiently dominant over the other terms. We refer to these as special 
regions. Outside these special regions, the terms remain bounded. Therefore, 
the ball resembles a piece of Swiss cheese, the holes (which may not be circular 
in general but in this case turn out to be nearly circular) being the special 
regions where movable singularities reside and the solid cheese being free of 
any singularity. 
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Starting at Zi in (the cheese-like region of) the ball, we continue the 
solution U along a ray until we encounter a point Z2 on the edge of a special 
region. Inside the region, we convert Pi to an integral equation by operating 
successively on the equation as though only the dominant terms were present. 

The dominant terms of 

u" = 6u^ + z, (2.32) 

are u" and 6m^. Integrating these dominant terms after multiplying by their 
integrating factor u', we get 

— = 2u^ + zu- udz + k, (2.33) 

2 J Z2 

with _ 

k := udz + k, 

J Z\ 

where the constant k (kept fixed below) is determined explicitly by the initial 
conditions. 

Since u is large, u' does not vanish (according to Eqn (|2.33|) ) and, there- 
fore, there is a path of steepest ascent starting at 22- We will use this idea 
to find a first point in the special region where u becomes infinite. 

Let d be an upper bound on the length of the path of integration from Z2 
to z and assume A > is given such that 

> AB, > AnB, A^ > 4d, A^ > 4\k\. 

Assume that \u\ > A a.t Z2. Then Eqn(^]3^) gives u'{z2) 7^= 0. Taking the 
path of integration to be the path of steepest ascent, we can show that 



\u 



'I > v^|n|3/2, 



and that the distance to a point where |m| becomes infinite is 

d < V2A-'/^ 



(See page 193 of for details.) So there is a first singularity encountered 
on this path which we will call zq. 

Now integrating the dominant terms once more (by dividing by 2u^, tak- 
ing the square root and integrating from zq) we get 

u= ( r h + — fk + zu- [\dz)] dz] . (2.34) 



2m3 



Z2 
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Substituting a function of the form 



where f{z) is analytic at zq into the right side of equation (p.34|) returns a 



function of the same form. Notice that, therefore, no logarithmic terms can 
arise. 

It is worth noting that the iteration of the integral equation( f^.34D gives 



the same expansion that we would have obtained by the Painleve test. In 
particular, it generates the appropriate formal solution without any assump- 
tions of its form and it point outs precisely where logarithmic terms may arise 
without additional investigations. (For example, try iteration of the integral 
equation with the term z on the right side of equation ( p. 32| ) replaced by z^, 
i.e. with zu — J^^ u dz replaced by z'^u — 2 J^^ zu dz in Eqn( p.34]) .) 



The remainder of the proof is a demonstration that the integral equation 
(|2.34| ) has a unique solution meromorphic in a disk centred at zq, that its 



radius is lower-bounded by a number that is independent of zq^ and that the 
solution is the same as the continued solution U . The uniformity of the lower 
bound is crucial for the proof. Uniformity excludes the possibility that the 
movable poles may accumulate to form movable essential singularities as in 
example (|2.31|) . 

Since the analytic continuation of U is accomplished along the union of 
segments of rays and circular arcs (skirting around the boundaries of succes- 
sive special regions encountered on such rays) and these together with the 
special regions cover the whole ball \z\ < B, we get a proof that the first 
Painleve transcendent is meromorphic throughout the ball. 

2.6 Rigorous Results for PDEs 

Sufficient results for the Painleve property of PDEs have been harder to 
achieve than ODEs. This is surprising because such results are lacking even 
for the most well known integrable PDEs. In this section, we describe some 
partial results towards this direction for the KdV equation. 

Definition 2.2 The WTC data for the KdV equation is the set {C,{t), U4{t), 
UQ{t)} of arbitrary functions describing this Painleve expansion. 



The following theorem proves that the series ( p. 301 ) converges for analytic 
WTC data. 
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Theorem 2.1 (Joshi and Petersen W^ ) Given an analytic manifold 
S := {{x,t) : X = ^{t)}, with ^(0) = 0, and two arbitrary analytic func- 
tions 

there exists in a neighbourhood of {0,0) a meromorphic solution of the KdV 
equation of the form 



w{x, t) = 777772 + ^(-^^ ^) 



-12 
(x^t)) 

where h is holomorphic. 



The next theorem provides us with a useful lower bound on the radius of 
convergence of this series. 



Theorem 2.2 (Joshi and Srinivasan j^^) Given WTG data ^{t), U4{t), 
Ueit) analytic in the ball -B2p+e(0) = {t G C : |t| < 2p + e}, let 

M= snp{l,\at)\,\u^{t)\,\u,{t)\}. 

\t\=2p 

The radius of convergence Rp = R of the power series (1-5) satisfies 

^ ^ min{l,p} 
- lOM 

This lower bound is used in to prove the well-posedness of the WTC 
Cauchy problem. That is, the locally meromorphic function function de- 
scribed by the convergent series ( ^.3U| ) has continuous dependence on the 
WTC data in the sup norm. 

To date there is no proof that the Korteweg-de Vries equation possesses 
the Painleve property. The main problem lies in a lack of methods for ob- 
taining the global analytic description of a locally defined solution in the 
space of several complex variables. However, some partial results have been 
obtained. 

The usual initial value problem for the KdV equation is given on the 
characteristic manifold t = 0. Well known symmetry reductions of the KdV 
equation (e.g. to a Painleve equation) suggest that a generic solution must 
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possess an infinite number of poles. WTC-type analysis shows that these 
can occur on noncharacteristic manifolds which intersect t = transversely. 
These results suggest that only very special solutions can be entire (i.e. have 
no singularities) on t = 0. 



Joshi and Petersen 371 showed that if the initial value 



u{x, 0) = uo{x) := a„a;", 

n=0 

is entire in x and, moreover, the coefficients real and nonnegative then 

there exists no solution holomorphic in any neighbourhood of the origin in 
unless 

Uq{x) = Oq + aix. 

This result can be extended to the case of more general under a condition 
on the growth of the function Uo{x) as a; — *• oo. 



3 Nonlinear-Irregular-Singular Point Analy- 
sis 

The Painleve expansions cannot describe all possible singular behaviours of 
solutions of differential equations. In this section, we describe some exten- 
sions based on irregular- singular-point theory for linear equations. 

The Painleve expansions at their simplest are Laurent series with a lead- 
ing term, and can, therefore, fail to describe solutions that possess movable 
isolated essential singularities. Consider the ODE 



3u'u"' = 5(u"? - (u''^ 
which has the general solution 

u{z) = a exp — ^o) ""^^^1 • 

Clearly u has a branched movable essential singularity. As suggested in ^ 
the Painleve test can be extended to capture this behaviour by considering 
solutions that become exponentially large near zq. To do this we expand 

u{z) = a_,{z)e^iz) + ao(z) + ai(z)e-^(^) + a2(z)e-2^(^) + a^{z)e-'-''^'^ + ..., 



u" iu'Y 



u 
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where S and the generahzed power series and S grows faster than any 

logarithm as z approaches 2:0- 

In other words, generahzed expansions (those involving logarithms, pow- 
ers, exponentials and their compositions) are necessary if we are to describe 
all possible singularities. These are asymptotic expansions which may fail to 
converge. They are in fact asymptotic expansions. We show in this section 
that they can nevertheless yield analytic information about solutions. We 
illustrate this with two main examples. The first is the Chazy equation and 
the second a fourth-order equation studied by Bureau. 



3.1 The Chazy Equation 

In this subsection, we examine the Chazy equation 

y"' = 2yy"-3{yr. (3.35) 



This equation is exactly solvable through the transformation [|T^, 



where Ui and U2 are two independent solutions of the Hypergeometric equa- 
tion 

_ , d'^u /I 1 \ du u 
^(^-l)3TT+ U-^0 -37- 777 = 0- 



dt"^ V2 Q J dt 144 
Following the work of Halphen Chazy noted that the function z{t) 



maps the upper half t-plane punctured at 0, 1, and 00 to the interior of a 
circular triangle with angles 7r/2, 7r/3, and in the 2-plane (see, for example, 



Nehari [^). The analytic continuation of the solutions Ui and U2 through one 
of the intervals (0, 1), (1, cxd), or (— cxo, 0) corresponds to an inversion of the 
image triangle across one of its sides to a complementary triangle. Continuing 
this process indefinitely leads to a tessellation of either the interior or the 
exterior of a circle on the Riemann sphere This circle is a natural harrier 
in the sense that the solution can be analytically extended up to but not 
through it. 

We will see below that any solution of equation ( p.35| ) is single-valued 
everywhere it is defined. The general solution, however, possesses a movable 
natural barrier, i.e. a closed curve on the Riemann sphere whose location 
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depends on initial conditions and through which the solution cannot be an- 
alytically continued. 

Leading order analysis of equation (|3.35|) shows that near a pole, 



y 6(z - zo) \ or ?/ ~ A{z - zq) ^ 

where A is an arbitrary (but nonzero) constant. On calculating successive 
terms in this generalized series expansion we find only the exact solution 

y{^) = I ^ .2 - — (3-36) 

[Z - ZqY Z- Zq 

This solution has only two arbitrary constants, A and Zq, and clearly cannot 
describe all possible solutions of equation (p. 35 ) which is third order. That 
is, solutions of the form ( |3.36| ) fail to capture the generic behaviour of the 
full space of solutions. 

The absent degree of freedom may lie in a perturbation of this solution. 
Applying the usual procedure for locating resonances i.e. substituting the 
expression 



|/(.)~-3^ + --- + /5(z-.o 



.r-2 



into equation (|3.35|) and demanding that {3 be free, we find 

(r + l)(r + 2)(r + 3) = 0. 

I.e. we must have r = —1, r = —2 or r = —3. The case r = —1 corresponds 
to the fact that A is arbitrary in ( 3.36 ). The case r = —2 corresponds to the 
freedom in zq. The case r = —3, however, indicates something more. 

Since the usual Frobenius-type series fails to describe the general solution 
near a singular point, we turn to be a nonlinear analogue of irregular singular 
point theory. Arguing from analogy with the linear theory (see, for example. 
Bender and Orszag |P) we seek a solution of the form 

y^'^ = (7^ - + ^'-'^^ 

where exp S{z) is regarded as small in a region near zq (generically, Zq will 
be on the boundary of this region). 
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Since the Chazy equation is autonomous we can, without loss of gener- 
ahty, take Zq = 0. For simphcity we also take A = 1/2. Substituting the 
expansion ( |3.37| ) into equation (|3.35|) gives 



S'" + 3S'S" + S'^ 
1 12 



S" + S"]+6^-^]S' 



z'^ z J ^ ' \z^ z^ 

+6 (^^ - ^) + 2 (S" + S'") e^ - 'iS'W (3.38) 

To ensure that exp(S') is exponential rather than algebraic, we must assume 
that S" <C S"^, 5"" <^ 5"^. Using these assumptions along with exp^* -C 1, 
S" ^ 1, and 2 -C 1, equation ( |3.38|) gives 

Z^ z 

Integration yields 

/X 1 6 k 

y{z) ~ TTT - - + -e ' . 
Iz^ z 

where k is an arbitrary constant; k represents the third degree of freedom 
that was missing from the Laurent series ( P^^Bj ). Extending to higher orders 
in expS'(z), we obtain the double series 



1 6 
2? ~ ^ 

+ 4e-'/^ (1 + 0{z)) + ^e-2/- (1 + 0{z)) (3.39) 
z-^ 82;^ 



+0 ^ 



Z2 



It can be shown that this series is convergent in a half-plane, given here by 
3fJ(l/2) > 0. 

This asymptotic series is valid wherever 

|fcexp(-l/2)| < 1. (3.40) 

Suppose k is small. Put 
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where ^ > to see whether the half-plane of validity can be extended. Then 
the condition (|3.40|) becomes 



< log 




By completing squares (after multiplying out the denominator) this can be 
rewritten as 

(e - 6r + » 5^ 

where 

5 := > 0. 

21og|A;| 

So asymptotically the region of vahdity of the series (|3.39|) lies outside a circle 
of radius 6 centered at —6. This is the circular barrier present in the general 
solution of the Chazy equation. In summary, the exponential (or WKB-type) 
approach has led to a three parameter solution. Morover, this description is 
valid in a region bounded by a circular curve where it diverges. 



3.2 The Bureau Equation 

Bureau partially extended the classification work of Painleve and colleagues 
to fourth-order equations. However, there were cases whose Painleve prop- 
erty could not be determined within the class of techniques developed by 
Painleve's school. One of these was 

u"" = 3u"m - 4u'^ (3.41) 

which we will call the Bureau equation. In this subsection, we show that the 
general solution is multivalued around movable singularities by using expo- 
nential or WKB-type expansions based on irregular-singular-point theory. 

It has been pointed out by several authors that Eqn (|3.41| ) possesses two 
families of Painleve expansions 

u^^ (3.42) 

(where we have shifted z — zq to z hj using the equation's translation invari- 
ance) distinguished by 

z/ = 2, 02 = 60 

u = 3, 03 arbitrary 



30 



with resonances given by 



z" 

z/ = 2 r = -3,-2,-1,20 
z/ = 3 r = -1,0. 

The case z/ = 2 has a full set of resonances (even though two are negative 
resonances other than the universal one). However, the case u = 3 is defec- 
tive because its perturbation (in the class of Painleve expansions) yields no 
additional degrees of freedom to the two already present in 03 and zq. It is, 
in fact, given by the two-term expansion 

u- ^ + 

z'^ 

We concentrate on this defective expansion in the remainder of this sub- 
section. Since this expansion allows no perturbation in the class of conven- 
tional Painleve expansions (which are based on regular- singular-point the- 
ory), we turn to perturbations of the form based on irregular-singular-point 
theory. Consider 

a 60 . 
u = — + — + u, 
z-^ z^ 

where 

u = exp{S{z)), S" > -, 2; < 1. 

(The assumption on S' is to assure that exp(S') is exponential rather than 
algebraic.) Substituting this into Eqn (|3.41|) , we get 



^'4 ^ ^s"+mS" + ^S'-QS''S" 
z-^ \ z-^ 

+ (^ + l^S--AS'S" 
\ z° z^ 




z^ j z^ 



+ (35" - S")e^ (3.43) 

The condition that exp(S') not be algebraic also implies that 5"^ ^ S", 
S'^ > S'", and S"^ > S'-^^l So dividing EgndOSp by S'^, taking the square 
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root of the equation, and expanding we get 

S' = ^ + g + 0(.-'/^) (3.44) 

where we have used recursive substitution of the leading values of S' and S" 
to get the term of order 1/z. That is, we get 

S = + — log z + const + o 1 . 

Take one such solution, with say a = 1/3. Then the perturbed solution has 
expansion 

n = ^ + 5 + fc±z^^/^exp(-2/zV2)(i + (3.45) 

where k± is an arbitrary constant. Note that there are two exponentials 
here (due to the two branches of the square root of z) and, therefore, k± 
represents two degrees of freedom. In the following, we consider only one 
of these solutions by fixing a branch of the square root in the exponential, 
say to be the one that is positive real on the positive real semi-axis in the 
2;-plane. For short, we write /c+ = k. 

Now consider the domain (or sector) of validity of this solution. Note 
that the neglected terms in its expansion contain a series of powers of exp(S') 
due to the nonlinear terms in Eqn (|3.43|) . Therefore, for the expansion to be 
asymptotically valid, this exponential term must be bounded i.e. 



kz^^/^e^v{-'^/z'"^) 



< 1 



^ \k\ e^^[^[-2/z'l' + (31/4) Hz))j < 1 (3.46) 

We show below that the domain of validity given by this inequality contains 
a punctured disk (on a Riemann surface) whose angular width is larger than 
27r. 

Assume there is a branch cut along the negative semi-axis in the 2;-plane 
with arg(z) G (— 7r,7r]. Consider z^^^ in polar coordinates, i.e. z^^^ = re*^, 
where — 7r/2 < 9 < -k /2. The positive branch will then give real part 

^(-2/z^'^ + (31/4) \n{z)) = -2 ^ cos(^^) + ^l^L, 
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Let K ■.= \n\k\/2. To satisfy Eqn(^.4^), we must have 



for r -C 1, i.e. 



2 311nr 
- - cos(^) + + 2K + oil) < 



cos(^)<-^-- Kr(l + o{l)). (3.47) 



Since r is small (note Inr < 0), this can only be violated near 9 = ±7r/2. 
Fix r small. Expand 6 = 7t/2 + e. Then Eqn( |3.47| ) gives 

e<-^^^^-Kr(l + o{l))+0{e^). (3.48) 



In particular, e can be negative (so long as |e| < 1). A similar calculation 
can be made near — 7r/2. 

These results show that the asymptotic validity of the solution given by 
Eqn( p.45D can be extended to a domain which is a disk of angular width 



> 27r. The small angular overlap is given by a sector of angular width 2e 
where e is bounded by O(rlnr) according to Eqn( |3.48| ). 

Let Zs be a point in this overlapping wedge with small modulus. At 
such a point, we have two asymptotic representations of u, one given by a 
prior choice of branch of z^^ and the other given by analytic continuation 
across the branch cut. If the true solution is single- valued in this domain, 
the choice of two asymptotic representations violates uniqueness. Therefore, 
the true solution must itself be multivalued. In other words, the exponential 
expansion shows that Bureau's equation cannot have the Painleve property. 



4 Coalescence Limits 

In this section we examine asymptotic limits of integrable equations that 
preserve the Painleve property. In the case of ODEs, such limits form the 
basis of Painleve's a-test. They are useful in the identification of noninte- 
grable equations and may be useful for indentifying new integrable equations 
as limits of others. 



Painleve noted that under the transformation 

u = ey + 



a = 4e 



-15 
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Pji becomes 



y"{x) 



Qy"^ + X + [2y^ + xy] . (4.49) 



In the limit as e vanishes, equation ( [4.49| ) becomes Pj. We write the above 



hmiting process as Pn — > Pi. Painleve gave a series of such asymptotic 
hmits which are summarized in Figure 1. 



/// 



Pvi — ^ Pv Pii — ^ Pi 



IV 



Figure 1: Asymptotic hmits among the Painleve equations. 

Each of these asymptotic limits coalesce the singular w-values of the 
Painleve equation (see Table 2) i.e. they coalesce movable singularities. In 



2^ , Halburd and Joshi proved that in the Pn — > P/ limit, simple poles of 
opposite residue coalesce to form the double poles in solutions of P/. They 
also proved that all solutions of Pj can be obtained as limits of solutions of 
equation (|]?S|): 

Theorem 4.1 Choose xo,a,P G C. Let yi and y he maximally extended 
solutions of Pi and equation ( \4-4^) respectively, both satisfying the initial 
value problem given by 

y{xo) = a, y'{xo) = 13. 

Let Q C C be the domain of analyticity of yj. Given any compact K G 

> such that y is analytic in (x, e) for x E K and |e| < rx- Moreover, 
y ^ yi in the sup norm as e — > 0. 

The series of asymptotic (or coalescence) limits given in Figure 1 by no 
means represents a complete list of such limits. The singular w-values of Pjv 
are and oo, corresponding to zeros and poles of the solutions respectively 
(the solutions of P/y are meromorphic |5^, ^). The standard coalescence 



limit in which P/y becomes Pn merges poles and zeros. However, the general 
solution of P/y contains simple poles of oppositely signed residue which may 
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be able to merge pairwise. An asymptotic limit coalescing these simple poles 
to form double poles does exist [p^ . 

To see this, consider a transformation in which regions near infinity (where 
the poles are close to each other) are mapped to the finite plane in the limit 
e — * 0. It is necessary to rescale u to counter a cancellation of the oppositely 
signed poles. This leads to new variables x and w{x) given by 

u{z) = e^w{x), 
z = N + e%, 

where p, g > 0, e <^ 1 and 1 is to be found in terms of e. Then Pjv 

becomes 

w,, = ^ + -e^^P+i^w'^ + ANeP+^'iw^ + 8eP+^'^xw^ + 2(N^-a)e^'iw 
2w 2 ^ ^ 

w H . (4.50) 

w 

The only maximal dominant balance {i.e. a limiting state of the equation in 
which a maximal number of largest terms remains P]) occurs when 



q = p, and a =: A^^ + ae 



-2q 



where a is a constant. Then setting the largest terms Ne''^'^^ and A^e^'^ to 
unity gives A^ = e~^^. Without loss of generality redefine ^ e. Then we 
get 

A^ = e"^ and a = e"^ + ae"^. 
Equation ( [4.50|) then becomes 

Wxx = — 4w;^ + (2x - a)w + — + -e'^w^ + Se'^xw'^ + 2e^x'^w, 
2w w 2 

or, in the limit e ^ 



w^^ = ^4:W^ + 2x - a)w + -, 
2w w 



which is Equation (XXXIV) (see p. 340 of Ince ||2^) in the Painleve- Gambler 
classification of second-order differential equations (after a simple scaling and 
transformation of variables). 
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Coalescence limits exist among PDEs also. For example, it is straightfor- 



ward to derive the transformation 

r = t; 



3 



u{x, t) 



which maps the modified Korteweg-de Vries equation (mKdV) 



Ut 



0, 



to 



Ur - 6e^U^U^ + 6UU^ + %5 = 0, 

which becomes the usual KdV equation in the limit e — > 0. An alternative 
method for obtaining the above asymptotic limit is to use the Pu — > Pj limit. 
The mKdV equation is invariant under the scaling symmetry 

X^^u, t ^— > X^t, X I— > Ax. 



u 

Define the canonical variables 

X 1 



z 



w 



logt, u{x,t) = i3t)-'/'y{z,w). 



(3t)i/3' - 3 
In terms of these variables, the mKdV equation becomes 



Vzz - 2?/ - zy - a)z + yy 



0, 



Pii 



where we have included the constant a to emphasize its relation to Pu. 
Now apply the asymptotic transformation used in the Pu — > Pj limit, to 
determine how y and z transform, and transform w in such a way that it 
remains in the limiting form of the equation as e — > 0. In this way we arrive 
at an equation equivalent to the KdV equation, which has a reduction to Pj. 



In ||2^ it is shown that the system 



Ex 
Ex 
2Nt 

Pt 
Pt 



P, 

-ipE + pE), 

NE, 

NE, 



(4.51) 
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admits a reduction to the full P/// {Pm with all four constants S 0, a, 
/3, 7 arbitrary). We note that if E = E* and p = p*, where a star denotes 
complex conjugation, we recover the unpumped Maxwell-Bloch equations; 

E, = p, 2Nt + pE* + p*E = 0, pt = NE. 

Consider solutions of system ( |4.51| ) of the form 

E = r'^e{z)w, 

E = t-^6{z)w-\ 

N = n{z), 

p = r{z)w, 

p = f{z)w^'^, 

where z := \fxi and w := (x/t)'^', k constant. Then 

y{z) := — - 

zr{z) 

solves Piii with constants given by 

a = 2{re-fe + Ak), (3 = A{l + 2k), 7 = 4(n^ + rf), 5 = -4. 
Note that by rescaling y we can make 5 any nonzero number. 



Using the procedure outlined for mKdV KdV, it can be shown that 
the Pjii Pii limit induces an asymptotic limit in which the generalized 
unpumped Maxwell-Bloch system ( [4.51| ) becomes the dispersive water-wave 
equation (DWW) 

Uxxxx + '^UtUxx + 4:UxUxt + QuiUxx + Uft = 0, 



which is known to admit a reduction to P// (Ludlow and Clarkson ||4^). The 
Pji Pj limit then gives DWW — > KdV. Ludlow and Clarkson have 
shown that DWW also admits a symmetry reduction to the full Piv- The 
PiY Pjj limit then induces an asymptotic limit that maps DWW back to 
itself in a nontrivial way. The limit Pjy — P34, outlined above, induces an- 
other limit in which DWW is mapped to the KdV. All six Painleve equations 
are known to arise as reductions of the self-dual Yang-Mills (SDYM) equa- 



tions (Mason and Woddhouse |^). Asymptotic limits between the Painleve 
equations can be used to induce similar limits within the SDYM system 
(Halburd 



37 



5 Acknowledgements 



The work reported in this paper was partially supported by the Australian 
Research Council. The authors also gratefully acknowledge with thanks the 
efforts of the organizing committee, particularly Dr K. M. Tamizhmani, and 
CIMPA in arranging the winter school on Integrable Systems at Pondicherry. 



References 

[1] Ablowitz, M.J. and Clarkson, P.A. (1991) "Solitons, Nonlinear Evolution 
Equations and Inverse Scattering," Lond. Math. Soc. Lecture Note Series 
149, Cambridge University Press, Cambridge. 

[2] Ablowitz, M.J. and Segur, H (1981) "Solitons and the Inverse Scattering 
Transform," SIAM, Philadelphia. 

[3] Ablowitz, M.J., Ramani, A. and Scgur, H. (1978) Nonlinear Evolution 
Equations and Ordinary Differential Equations of Painleve Type, Lett. 
Nouvo Cim., 23, 333-338 

[4] Ablowitz, M.J., Ramani, A. and Segur, H. (1980) A Connection between 
Nonlinear Evolution Equations and Ordinary Differential Equations of 
P-type. I and II, J. Math. Phys., 21, 715-721 and 1006-1015 

[5] Ablowitz, M.J. and Segur, H. (1977) Exact Linearization of a Painleve 
Transcendent, Phys. Rev. Lett. 38, 1103-1106 

[6] Arnold, V.I. (1984) Mathematical Methods of Classical Physics, 
Springer- Verlag, New York 

[7] Abramowitz, M. and Stcgun, I. eds. (1972) Handbook of Mathematical 
Functions, Dover, New York. 

[8] Bender, CM. and Orszag, S.A. (1978) "Advanced Mathematical Meth- 
ods for Scientists and Engineers," McGraw-Hill, New York 

[9] Bountis, T., Segur, H. and Vivaldi, F., (1982) Integrable Hamiltonian 
systems and the Painleve property, Phys. Rev. 25A, 1257-1264. 



38 



[10] Boutroux, P. (1913) Recherches sur les transcendents de M. Painleve et 
r etude asymptotique des equations differentielles du second ordre, Ann. 
Ecole Norm. Super., 30, 255-375; 31, 99-159 

[11] Bureau, F.J. (1964) Differential equations with fixed critical points, An- 
nali di Matematica pura ed applicata LXVI, 1 - 116. 

[12] Chazy, J. (1910) Sur les equations differentielles dont I'integrale generale 
possede une coupure essentielle mobile, C.R. Acad. Sc. Paris, 150, 456- 
458 

[13] Chazy, J. (1911) Sur les equxitions differentielles du troisieme et d'ordre 
superieur dont Vintegrale generale a ses points critiques fixes. Acta 
Math., 34, 317-385 

[14] Coddington, E. A. and Levinson, N. (1955) Theory of Ordinary Differ- 
ential Equations, McGraw-Hill, New York. 

[15] Clarkson, P. A., and Cosgrove, C. M., (1987) Painleve analysis of the 
non-linear Schro dinger family of equations J. Phys. A: Math. Gen. 20, 
2003-2024. 

[16] Conte, R., Fordy, A., and Pickering, A., (1993) A perturbative Painleve 
approach to nonlinear differential equations, Physica 69D 33-58. 

[17] Cosgrove, C. M., (1993) Painleve classification of all semilinear PDEs of 

the second order, I Hyperbolic equations in two independent variables, II 
Parabolic and higher dimensional equations. Stud. Appl. Math.89 1-61, 
95-151. 

[18] Fordy, A and Pickering, A. (1991) Analysing Negative Resonances in the 
Painleve Test, Phys. Lett. 160A, 347-354 

[19] Fuchs, R (1907) Math. Annalen 63 301-321 

[20] Gambier B. (1910) Sur les equations differentielles du second ordre et du 
premier degre dont I'integrale generale est a points critiques fixes. Acta 
Math., 33, 1-55 

[21] Gardner, C.S., Greene, J.M., Kruskal, M.D. and Miura, R.M. (1967) 
Method for Solving the Korteweg-de Vries Equation, Phys. Rev. Lett., 
19, 1095-1097 



39 



[22] Grimshaw, R. (1979) Slowly Varying Solitary Waves. I. Korteweg- 
de Vries Equation, Proc. R. Soc. Lond. 368A, 359-375 

[23] Halburd, R. (1996) Integrable Systems; their singularity structure and 
coalescence limits, PhD thesis, University of New South Wales 

[24] Halburd, R. and Joshi, N. (1996) The Coalescence Limit of the Second 
Painleve Equation, (to appear in Stud. Appl. Math.) 

[25] Halburd, R. and Joshi, N. (1996) Coalescence Limits of Integrable Partial 
Differential Equations, in preparation. 

[26] Halphen, G. (1881) Sur un systeme d'equations dijferentielles, 
C.R. Acad. Sc. Paris, 92, 1101-1103 and 1404-1407 [(Euvres (1918), 
Paris, 2 475-477 and 478-481] 

[27] Ince, E.L. (1956) "Ordinary Differential Equations," Dover, New York 

[28] Hlavaty, L., The Painleve analysis of damped KdV equation, J. Phys. 
Soc. Japan 55 (1987) 1405-1406. 

[29] Hlavaty, L. , (1990) Painleve classification of PDEs, in Proc. of NATO 
Advanced Research Workshop, P. Winternitz and D. Levi (eds), Ste 
Adele, Quebec Sept 2-7, 1990, (Plenum Publishing Co., 1991). 

[30] Joshi, N., (1987) Painleve property of general variable-coefficient ver- 
sions of the Korteweg-deVries and nonlinear Schrddinger equations, 
Phys. Letts. 125A, 456-460. 

[31] Joshi, N. and Kruskal, M. (1992) A Direct Proof that Solutions of the 
First Painleve Equation have no Movable Singularities except Poles, in 
"Nonlinear Evolution Equations and Dynamical Systems" (Baia Verde, 
1991) 310-317, Boiti, M., Martina, L., and PemineUi, F. (eds) World 
Sci. Pubhshing, River Edge, NJ 

[32] Joshi, N. and Kruskal, M.D. (1993) A New Coalescence of Movable Sin- 
gularities in the Fourth Painleve Equation, University of New South 
Wales Apphed Mathematics Report AM93/11 

[33] Joshi, N. and Kruskal, M.D. (1993) A Local Asymptotic Method of See- 
ing the Natural Barrier of the Solutions of the Chazy Equation in "Ap- 
plications of Analytic and Geometric Methods to Nonlinear Differential 



40 



Equations (Exeter, 1992)," 331-339, NATO Adv. Sci. Inst. Ser. C Math. 
Phys. Sci., 413, Kluwer Acad. Publ., Dordrecht 

[34] Joshi, N. and Kruskal, M.D. (1994) A Direct Proof that Solutions of 
the Six Painleve Equations have no Movable Singularities Except Poles, 
Stud. Appl. Math., 93, 187-207 

[35] Joshi, N. and Petersen, J.A. (1996) Complex blow-up in Burgers' equa- 
tion: an iterative approach, Bull. Austral. Math. Soc. 54, 353-362 

[36] Joshi, N. and Petersen, J.A. (1994) A Method for Proving the Con- 
vergence of the Painleve Expansions of Partial Differential Equations, 
Nonhnearity, 7, 595-602 

[37] Joshi, N. and Petersen, J.A. (1996) Nonexistence Results Consistent with 
the Global Painleve Property of the Korteweg-de Vries Equation, (to 
appear) 

[38] Joshi, N. and Srinivasan, G.K. (1997) The Radius of Convergence and 
the Well-Posedness of the Painleve Expansions of the Korteweg-de Vries 
Equation, Nonhnearity, to appear. 

[39] Kowalevski, S. (1889a) Sur le Probleme de la Rotation d'un Corps Solid 
Autour d'un Point Fixe Acta Math., 12, 177-232 

[40] Kichenassamy, S., and Littman, W., (1993) Blow-up surfaces for non- 
linear wave equations I, II Comm. PDEs 18 431-452 and 1869-1899. 

[41] Kichenassamy, S., and Srinivasan, G., (1995) The structure of WTC 
expansions and applications J. Phys. A 28 1977-20004. 

[42] Kowalevski, S. (1889b) Sur une Propriete d'un Systeme d'Equations 
Differentielles qui Definit la Rotation d'un Corps Solide Autour d'un 
Point Fixe, Acta Math., 14, 81-93 

[43] Kruskal, M.D. (1963) Asymptotology, in "Mathematical Models in Physi- 
cal Sciences" Dobrot, S. ed, Prentice-Hall, Englewood Cliffs, New Jersey, 
17-48. 



41 



[44] Kruskal, M.D. (1992) Flexibility in Applying the Painleve Test, in 
"Painleve Transcendents - their Asymptotics and Physical Apphca- 
tions" Levi, D. and Winternitz, P. eds NATO ASI Series B: Physics 
vol 278, Plenum Press, New York 

[45] Kruskal, M.D. and Clarkson, P.A. (1992) The Painleve-Kowalevski and 
Poly-Painleve Tests for Integrability , Stud. Appl. Math., 86, 87-165 

[46] Kruskal, M. D., and N. Joshi, N. (1991) Soliton theory, Painleve property 
and integrability, in Chaos and Order, Miniconference on. Proceedings 
of the Centre for Mathematical Analysis, Australian National Univer- 
sity, N. Joshi and R. L. Dewar (eds) Feb. 1-3, 1990, (World Scientific 
Publishing Co., Singapore) 82-96. 

[47] Ludlow, D.K. and Clarkson, P.A. (1993) Symmetry Reductions and Ex- 
act Solutions for a Generalised Boussinesq Equation in "Applications of 
Analytic and Geometric Methods to Nonlinear Differential Equations" 
P.A. Clarkson (Ed.) NATO ASI Series C vol 413, Kluwer, Dordrecht, 
415-430 

[48] Mason, L.J., and Woodhouse, N.M.J. (1993) Self-Duality and the 
Painleve Transcendents Nonlinearity, 6, 569-581 

[49] McLeod, J.B. and Olver, P.J. (1983) The Connection Between Partial 
Differential Equations Soluble by Inverse Scattering and Ordinary Dif- 
ferential Equations of Painleve Type, SIAM J. Math. Anal., 14 , 488-506 

[50] Nehari, Z. (1952) "Conformal Mapping," McGraw-Hill, New York 

[51] Newell, A. C, Tabor, M., and Zheng, Y., (1987), A unified approach to 
Painleve expansions Physica 29D, 1-68. 

[52] Olver, F.W.J. , (1992), "Asymptotics and Special Functions," Academic 
Press, London. 

[53] Painleve, P. (1888) Sur les equations differentielles du premier ordre, 
C.R. Acad. Sc. Paris, 107, 221-224, 320-323, 724-726 

[54] Painleve, P. (1900) Memoire sur les equations differentielles dont 
I'integrale generale est uniforme. Bull. Soc. Math. France, 28, 201-261 



42 



[55] Painleve, P. (1902) Sur les equations differentielles du second ordre et 
d'ordre superieur dont I'integrale generale est uniforme, Acta Math., 25, 
1-85 

[56] Painlcvc, P. (1906) Sur les equations differentielles du second ordre d 
points critiques fixes, C. R. Acad. Sc. Paris, 143, 1111-1117 

[57] Ramani, A. , Dorizzi, B. and Grammaticos, B., (1984) Integrability and 
the Painleve property for low- dimensional systems, J. Math. Phys. 25, 
878-883. 

[58] Ramani, A., Grammaticos, B., and Bountis, T., (1989) The Painleve 
property and singularity analysis of integrable and non-integrable sys- 
tems, Phys. Rep. 180, 159-245. 

[59] Weiss, J., Tabor, M. and Carnevale, G. (1983) The Painleve Property 
for Partial Differential Equations, J. Math. Phys., 24, 522-526 

[60] Yoshida, H. (1983), Necessary condition for the existence of algebraic 
first integrals: I, II, Celes. Mech. 31, 363-379; 381-399. 

[61] Ziglin, S. (1983), Branching of solutions and nonexistence of first inte- 
grals in Hamiltonian mechanics; I, II, Func. Anal. Appl. 16 181-189; 
17, 6-17. 

[62] Ziglin, S. (1982), Self-intersection of the complex separatrices and the 
nonexistence of the integrals in the Hamiltonian systems with one-and- 
half degrees of freedom, J. Appl. Math. Mech. 45 411-413. 



43 



